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Abstract 

A nonequilibrium radiative heating prediction method has been used to evaluate several 
energy exchange models used in nonequilibrium computational fluid dynamics methods. The 
radiative heating measurements from the FIRE II flight experiment supply an experimental 
benchmark against which different formulations for these exchange models can be judged. 
The models which predict the lowest radiative heating are found to give the best agreement 
with the flight data. Examination of the spectral distribution of radiation indicates that 
despite close agreement of the total radiation, many of the models examined predict excessive 
molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may 
lead to a correction for this problem. 


Introduction 

Many proposed missions envision the use of aerobraking to return payloads from geo- 
stationary orbit, the Moon, and Mars. To design vehicles for this purpose, knowledge of the 
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total heating environment is required. The radiative heat load on a vehicle moving through 
the atmosphere increases as the velocity and size increase, so that many of these missions 
will experience significant radiative heating. Accurate prediction of radiative as well as con- 
vective heating therefore becomes important to designers and mission planners. Since these 
vehicles will spend a large portion of their atmospheric trajectory at high altitude, where low 
freestream density can lead to significant thermal and chemical nonequilibrium effects, the 
tools used to predict the heating environment must account for nonequilibrium. Radiative 
heating predictions therefore depend on flowhelds obtained from nonequilibrium computa- 
tional fluid dynamics (CFD) codes, like LAURA (the Langley Aerothermodynamic Upwind 
Relaxation Algorithm 1 ). These flowheld solvers contain energy exchange models that have 
yet to be validated. While perturbations in these models may have little effect on convective 
heating, their effect on the radiative heating level can be large. Experimental data on this 
question are limited, consisting almost exclusively of the FIRE flights of the mid-1960’s. 2 

The objective of this study is to examine a number of the energy exchange models 
implemented in LAURA for their influence on stagnation point nonequilibrium radiative 
heating. By comparison with the FIRE II flight radiometer data, some insight into the 
validity of the various forms proposed for these exchange models may be obtained. 

LAURA is a nonequilibrium CFD code that uses Park’s two-temperature approximation, 3 
in which the heavy-particle-translational/rotational energy mode (hereafter referred to as the 
translational mode), and the vibrational/electronic/electron-translational energy mode (re- 
ferred to as the vibrational mode 4 ) are considered separately. LAURA contains a number of 
semi-empirical models governing the exchange of energy between these modes. 

The radiation model used in this study is the Langley Optimized Radiative Nonequi- 
librium (LORAN) code, which is a recent development. 5 It treats the radiation spectrum in 
moderate detail, accounting for all mechanisms of importance. In this paper, only the sub- 
range of the visible and infrared spectrum where the FIRE radiometer windows are transpar- 
ent is included. Gas absorption is negligible in this spectral region so the radiative transport 
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is calculated assuming an optically thin shock layer. The tangent slab approximation is also 
invoked since the FIRE forebody is blunt. 

The FIRE flight experiment was conducted during the 1960’s to determine the heat 
load on an Apollo-type vehicle entering Earth’s atmosphere at 11.4 km/sec. 2,6 Of the two 
flights, the first experienced telemetry problems that made data reduction and analysis 
difficult as well as control problems during the second half of the entry that resulted in 
substandard data. These problems were corrected before the second flight, which provided 
good measurements of the total and spectral radiation at the stagnation point during the 40 
second entry period. 7,8 The second flight will be emphasized here. 

Solutions are obtained with the LAURA code at three points in the FIRE II trajectory 
that span the range from extreme nonequilibrium early in the flight to a near-equilibrium 
point approaching peak heating. Each point is solved using the different energy exchange 
models under consideration. The flowheld conditions are then input to the LORAN code 
to obtain predictions of the radiative heating. The results are compared with the flight 
radiometer data to evaluate the validity of the various energy exchange models. 

Theory 

The LAURA code contains a number of semi-empirical energy exchange models that 
govern the transfer of energy between the translational and vibrational modes under nonequi- 
librium conditions. Several of these models continue to be the subject of debate 4 as there has 
been no definitive determination of their validity because of a lack of experimental or flight 
data. This study will examine the effect of models for the rate controlling temperature for 
dissociation, vibrational relaxation cross section, and energy exchange during dissociation. 
The effect of chemical kinetics is briefly examined as well. Energy exchange during electron 
impact ionization was considered, but since it has little effect under the low ionization levels 
of the FIRE flights, it was not pursued. 
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Dissociation Temperature 

The dissociation temperature, T d} is the rate-controlling temperature for dissociation 
reactions. Several formulations have been proposed for this temperature, notably by Park 
and his co-workers. 3,9 These models are empirical and consist of a geometric weighting of 
the translational temperature T t and the vibrational temperature T v as follows: 


T d 


rpmrpr 

1 v 


( 1 ) 


The choice of powers m and n adjusts the weight given to each temperature. The recent 
trend has been to increase the weighting of the translational temperature T t because of indi- 
cations that heavy particle collisions are more important than other mechanisms in causing 
dissociation. This study considers three different models with (m,n) equal to (.5, .5), (.7, .3), 
and (1,0). The first set is the original Park two-temperature proposal, the second the more 
recent model, and the third is included to study the extreme case of dissociation controlled 
by the heavy particle temperature alone. 

For convenience in identifying the various cases in this paper, a run identifier is used. 
The dissociation temperature models will be denoted by 55, 7 5, and 10 } respectively, in the 
third and fourth digits of this identifier. 


Vibrational- Translational Energy Exchange Cross Section 

The equilibration of the vibrational and translational energy modes is modelled using 
a relaxation time, t„, for each species. Millikan and White 10 proposed a semi-empirical 
formulation for T^f w in the range of 300 to 8000 K. Park 11 suggests an additional collision 
limiting correction, , for high temperatures where the Millikan- White correlation predicts 
excessively fast relaxation. The total relaxation time is then 


where Park’s contribution is 


Tv = T„ + T 


t„ p = (a„c,n s ) 1 
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In this expression cj is the average velocity and n s the number density of molecular species 
s, and (T v is an effective cross section for vibrational relaxation. The latter quantity has been 
the subject of some debate, with three different values receiving support in various efforts 
to match the limited experimental data that has bearing on the question: 


(50, 000/Tj) 2 m 2 

( 4 ) 

: 10“ 2 ° to 2 

(5) 

: 10“ 21 to 2 

(6) 


Equation 5 is the original proposal, which assumes this cross section to be one tenth of the 
elastic cross section. Equation 6 was suggested when Eq. 5 seemed still too high. Equation 4 
is a modification of Eq. 6 that reduces the contribution of r ^ at low temperatures. All 
three models are considered in this study, and are identified by the notation vl, v2 } and v3, 
respectively. 

A further modification to this energy exchange process has been proposed by Park. 3 
He argues that vibrational relaxation exhibits a diffusion-like behavior at high temperatures 
which requires a correction to the relaxation time. However, this modification requires 
evaluation of the post-shock levels of T t and T v . Interpretation of these quantities in a shock 
capturing solution is somewhat ambiguous. Improper definition can lead to instabilities in 
some circumstances, consequently no attempt was made to include this modification. 

Energy Exchange in Dissociation 

When a diatomic molecule dissociates, the vibrational energy it contained is consumed 
by the higher ground states (heats of formation) of the constitutive atoms. The energy thus 
removed must be accounted for in the vibrational energy equation. The amount of energy 
lost is commonly assumed to be the average vibrational energy at the local conditions, ejh 
Because dissociation from a higher vibrational state may be more probable (the concept of 
preferential dissociation), several other formulations have been proposed. 4 
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One model assumes that the vibrational energy loss in dissociation, A E v , is some fraction 
of the dissociation energy of a molecule measured from its ground state, D: 

A E v = c x D (0 < a < 1) (7) 

Sharma et al. 9 have proposed using c x = 0.3. This model was considered, but since it predicts 
the non-physical result of a negative T v in the shock region, it was not pursued. 

Park 12 has proposed that 

A E v = D — kT t (8) 

where k is Boltzmann’s constant, so that the energy lost is the dissociation energy minus the 
average translational energy. This model arises from the assumption that dissociation only 
occurs from levels that are within the average collisional energy of the dissociation threshold 
(similar to the mechanism of bound-free radiation). Again, this A E v is too large and results 
in an unrealistically small T v behind the shock. 

A third proposed model assumes that dissociation occurs from some vibrational state(s) 
above the average. This is expressed as: 


A E v = c 2 e v (c 2 > 1.) (9) 

The best value of c 2 to use in this empirical model is unknown, and may depend on e7 (or 
T v ). In this work c 2 = 2 has been selected to provide an initial assessment of the model. The 
suffix 2evs is used in the run identifier for these cases. 

FIRE II Flight Experiment 

The FIRE flight vehicle was an Apollo-like geometry with a layered beryllium heatshield. 
(The second layer had a geometry identical to Apollo.) Each of the three layers was used 
up to a temperature limit and jettisoned, resulting in three periods of prime data during 
the entry. The first heatshield had a nose radius of 0.935 m and a diameter of 0.672 m. All 
the cases considered here are for this first heatshield, which was jettisoned near the point of 
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Table 1: Selected FIRE II Trajectory Conditions 


Case 

Time 

Altitude 

Velocity 

Density 

T 

- L W 


(sec) 

(km) 

(km/sec) 

(kg/m 3 ) 

(K) 

A 

1631 

84.59 

11.37 

9.15e-6 

460 

B 

1634 

76.41 

11.36 

3.72e-5 

615 

C 

1637.5 

67.04 

11.25 

1.47e-4 

1030 


peak heating. For this study it has been modelled as a sphere with an effective nose radius 
of 0.747 m. A fully catalytic wall boundary condition was used. The cases selected from 
the FIRE II flight 6 are listed in Table 1. These conditions cover the range from extreme 
nonequilibrium to near-equilibrium 7 and were chosen to ensure that the energy exchange 
models are assessed in widely disparate flow regimes. The temperature profiles predicted 
by LAURA along the stagnation line for Case C are presented in Fig. 1, using the baseline 
set of energy exchange models denoted vl73. The temperature profiles for Cases A and B 
are shown with the label vl73 in Figs. 3 and 7. In these figures, y is the distance from the 
wall. Case A is in extreme thermal nonequilibrium, while Case C is in thermal equilibrium 
through about half the shock layer. 

The spectral response of the radiometer windows is reported by Dingeldein. 13 It is flat 
between 0.23 and 2 /mi in wavelength, falling off sharply below 0.2 /mi and more gradually 
above 2 /mi. The design goal is reported as a flat response from 0.2 to 6 /mi. The various 
reports on FIRE quote actual spectral ranges starting between 0.2 and 0.23 /mi, and ending 
between 4 and 6 /mi. Taking these variations into account, the spectral range of 0.23 to 
4 /mi (0.31 to 5.4 eV) is considered to be closest to the actual window transmission range, 
and is used here. Though little radiation occurs between 4 and 6 /mi at these conditions, the 
reduced spectral range is used to provide more exact comparisons. FIRE also had scanning 
spectral radiometers that were designed to cover the range of 0.2 to 0.6 /mi (2.1 to 4.1 eV). 
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Mechanical problems during the flight limited forward scans to 0.3-0.5575 /mi, and backward 
scans to 0.6090-0.3 /mi. 8 The spectral resolution is quoted as 0.004 /mi, with a root-sum- 
square uncertainty in the measured spectra of ± 23 percent. This resolution is insufficient 
to resolve any atomic lines. 

The radiation levels obtained from the flight data are given in Tables 3-5. These results 
were taken from Cauchon 7 and converted from intensities to fluxes assuming a transparent 
plane-parallel shock to be consistent with the present calculations. Also shown is the root- 
sum-square uncertainty estimated in a post-flight error analysis. It should be noted that the 
uncertainty in the radiation level for Case A is probably higher, as the radiometers were at 
the lower limit of their sensitivity range for that case. The measurements show considerable 
scatter during the early parts of the trajectory, with a variation of about a factor of three 
between high and low readings for the first several seconds. The scatter decreases as the 
level of radiation increases. It is noteworthy that the FIRE I data acquired in this density 
range for a slightly higher velocity (11.57 vs 11.37 km/sec) were higher by a factor of three 
than the FIRE II results for both the total and spectral radiometers. The two flights thus 
suggest a very sensitive dependence of radiative heating on freestream velocity under these 
flow conditions. At the conditions of Case B, the flight-to-flight variation is about a factor 
of 1.6. For Case C, it is about 1.3. 

Results and Discussion 

All the cases in this paper were run with 64 grid cells normal to the body in the 
LAURA solution. This is not sufficient to resolve the shock, especially in Cases B and 
C. One grid-resolved solution was obtained to quantify the effect, and indicated that the 
radiation prediction varies only about 15 percent. This is not enough to change any of the 
conclusions drawn from this study. 
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Table 2: Validation of Radiation Model 


Case 

Time 

(sec) 

Flux, W/cm 2 (0.31-6.2 eV) 

LORAN 

NEQAIR 

(Emission) 

Deviation 

(percent) 

A 

1631.0 

2.76 

2.47 

-11.7 

B 

1634.0 

22.2 

21.8 

-1.8 

C 

1637.5 

91.5 

83.0 

-10.2 


Radiation Model 

The initial validation of the LORAN nonequilibrium radiation model is reported in 
Ref. 5. To illustrate its accuracy for the FIRE II flow conditions, Table 2 compares it to the 
NEQAIR code of Park. 14 NEQAIR is a detailed line-by-line nonequilibrium radiation model. 
The NEQAIR results reported here were obtained using its “emission only” mode, which 
sums the total emission from each radiating line without considering the line shape. This 
mode is fast and avoids problems of spectral resolution. The agreement between the two 
methods is within 12 percent and is considered excellent. (The slight differences between 
this table and Tables 3-5 are caused by changes made to the LAURA code between the two 
sets of computations and by the use of different spectral limits.) 

Energy Exchange Models 

The radiation predictions for the various energy exchange models are summarized and 
compared with the flight data in Tables 3 to 5. The predictions are broken out into atomic 
continuum, molecular band, and atomic line contributions. Model variations are denoted by 
the run identifiers defined above in the discussion of the various energy exchange models. 
For example, the run identifier vl73 implies the use of Eq. 4 for the vibrational-translational 
energy exchange cross section, Td = TfT^ 3 7 and the use of the average vibrational energy 
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in the dissociation energy loss. Though the energy exchange processes examined are all 
interrelated, only one model is varied at a time from the vl73 model baseline. The effect of 
each energy exchange model and its comparison to the flight data are discussed in detail in 
the sections that follow. 

Dissociation Temperature 

As shown in the stagnation line temperature profiles, the translational temperature 
T t is much higher than T v near the shock. Increasing the weighting of the dissociation 
temperature, T^, on T t in Eq. 1 therefore results in faster dissociation of the molecular species. 
This dissociation removes both translational and vibrational energy, decreasing T t and T v . 
The increased dissociation is shown in Fig. 2 which compares the N 2 and 0 2 profiles along 
the stagnation line for Case A for the two extreme Td models. The effect on the temperature 
profiles is shown in Fig. 3 also for Case A. The decreased temperature results in increased 
density and a smaller shock standoff distance. For the three FIRE II cases, placing increasing 
weight on the translational temperature T t in the definition of the dissociation temperature 
Td decreases the total radiation. Part of the decrease occurs in the molecular bands, which 
depend directly on the molecular concentrations. The rest of the decrease results from the 
lower T v and decreased standoff distance affecting all three radiating mechanisms. Figure 4 
shows the emission profiles for the three Td models for Case A. The sensitivity to the Td 
model decreases at the later flight times of Cases B and C. 

Tables 3-5 reveal that in all three cases the predictions are closer to the flight data when 
T t is weighted more in the definition of Td } confirming other recent work on this model. For 
Cases A and B, agreement within about a factor of three is found. For Case C, the variations 
caused by this model are within the uncertainty of the flight data. 
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Vibrational-Translational Energy Exchange Cross Section 


Figures 5-6 present the effect of the different vibrational-translational energy exchange 
cross sections in Eqs. 4-6 on the temperature and radiative emission profiles. For a smaller 
<r.y, the relaxation time t v increases. Increasing t v delays the relaxation of translational to 
vibrational energy. This increases T t while reducing T„, as illustrated in Fig. 5 for Case A. 
Although the standoff distance increases with increasing T t , the reduction in radiation due to 
the lower T v in this model more than compensates. The lower T v reduces the radiation from 
all three mechanisms. Figure 6 shows the emission profiles for Case A. The lowest radiation 
for all three FIRE cases is predicted with Eq. 6, and the highest with Eq. 5. Comparing 
the predicted radiation for the vl73, v273, and v373 models in Tables 3-5 shows that the 
influence of o v is largest in Case A. 

Again all three cases are closest to FIRE II when the radiation prediction is lowest, with 
Eq. 6, though only Case A exhibits much sensitivity or discrimination among the models. The 
predictions in Case B remain high, while the effect on Case C is within the data uncertainty. 

Energy Exchange in Dissociation 

Figure 7 compares the temperature profiles for Eq. 9 with c 2 = 1 (the baseline model) 
and c 2 = 2 for Case B. The increased A E v in the 2evs model (c 2 = 2) reduces T v in the 
region behind the shock where dissociation occurs, and increases it deeper in the shock layer 
where recombination begins. It also results in a smaller shock standoff distance. Figure 8 
presents the radiative emission for Case B and confirms that the lower T v near the shock 
results in reduced radiation, while the higher T v in the layer increases it. This is observed in 
the tables as a decrease in the molecular band radiation and a slight increase in the atomic 
contributions, resulting in a net decrease in the total radiation. (The lower peak in the 
dashed curve of Fig. 8 is the atomic radiation peak.) 

For Cases A and B, the decreased radiation resulting from the 2evs model brings the 
prediction closer to the flight measurement, but it remains high. For Case C the effect is 
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within the data uncertainty, but this is the only model where the prediction is lower than 
the flight measurement. Table 5 also includes a solution run with the dissociation energy 
exchange model of Eq. 7 using c\ = 0.3 (denoted 3dis) for comparison. A flag was introduced 
in LAURA to suppress negative temperatures. No significant variation is observed that would 
encourage the use of this model. 

As with the two previous energy exchange models, the sensitivity to this model decreases 
as the flowheld tends toward equilibrium. The decreasing sensitivity to all these models near 
equilibrium is expected, since they govern the exchange between energy modes in the two- 
temperature thermal nonequilibrium model. When dealing with flowhelds near equilibrium, 
the choice of formulation for these models is not as important. 

Spectral Radiation Comparison 

To obtain further information on the validity of the energy exchange models, the FIRE 
spectral radiometer measurements (vs. wavelength, A, in /mi) are compared with those ex- 
change models whose predictions are closest to the flight measurements. Only continuum 
radiation is shown for the predicted spectra, since the resolution of the flight spectral ra- 
diometer is not sufficient to record any line radiation. In evaluating the flight spectra pre- 
sented in this paper, it should be noted that only a few are available in the literature. Those 
that have been published are the ones that were judged to contain the least amount of noise. 
This criterion had the effect of selecting spectra that correspond to the lower range of flight 
radiation measurements. It is not known whether the unpublished data which correspond to 
the higher range of measurements had increased intensity levels or broader spectra or both. 

Case A (t=1631 sec) 

The spectral distribution of radiation has been examined for Case A for the v313 
solution. Figure 9(a) shows the prediction, while Fig. 9(b) is a flight spectral radiometer 
scan at 1631.3 sec reproduced from Cauchon. 8 While the prediction is much higher than the 
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flight scan, they are qualitatively similar. Both show radiation from the iV^ first negative 
band, whose two major band heads can be distinguished. (The spike in the flight spectrum 
between 0.5 and 0.6 /mi has been attributed to a spurious signal. 8 ) Integrating the LORAN 
predictions over this limited spectral range results in a value that is high but within the data 
scatter at this flight condition. 

Case B (t=1634 sec) 

Figure 10 compares the 2evs result for Case B to a flight spectrum at 1634.43 sec. 8 
The qualitative agreement is comparable to that shown for Case A in Fig. 9. Integrating 
the predicted radiation over the range of the spectral radiometer gives a result that again, 
while high, is within the scatter of the flight data given in Fig. 13 of Cauchon. 7 Also shown 
on Fig. 10 is a spectrum from a solution using the Dunn & Kang chemical kinetics model. 15 
The radiation prediction for this case is included in Table 4, where it is denoted dk73. 
This prediction is much closer to the FIRE II measurements than any of those obtained for 
Case B with the Park kinetics. 3 (The same cannot be said in Case A (see Table 3).) The 
differences between these two chemical kinetics models are complex and should be studied 
further. It is interesting to note, however, that the distribution and relative magnitudes of 
the several band heads in this figure agree quite well with the flight data, and the Dunn & 
Kang prediction exhibits quantitative as well as qualitative agreement with the measured 
spectrum. 

Case C (t=1637.5 sec) 

Figure 11 compares the predicted v373 spectrum for Case C to the nearest available 
flight spectrum, at 1636.43 sec. By this point in the trajectory, the radiometer windows are 
approaching their melting point. A detailed post-flight analysis determined, however, that 
the windows cause less than a 10 percent change in the measured radiation. 7 The predicted 
spectrum is excessive. Because this case is close to thermal equilibrium, solutions were also 


13 



generated assuming a single temperature, for both the Park and Dunn & Kang chemical 
kinetics. These allow an assessment of the extent and impact of thermal nonequilibrium, 
and are included in Table 5 ( vleq and dkeq, respectively). The vleq result is shown on 
Fig. 11 as the dashed line. The agreement between this spectrum and the flight spectrum is 
remarkable. (Recall that the spectral radiometer did not measure below 0.3 /mr.) 

Cases A, B, and C 

Referring again to the temperature profiles for the three cases, it is clear that the 
thermal relaxation is in two stages with different time constants. Right behind the shock, 
T t decreases rapidly as a result of dissociation. The decrease in T t then slows indicating 
the involvement of other slower reactions, such as ionization. This two-stage equilibration is 
much less apparent for the Dunn & Kang chemical kinetics model (see Fig. 12), and is not 
an issue in the equilibrium solutions. The improved agreement noted for Case B with the 
Dunn & Kang solution and for Case C with the equilibrium solution suggests, therefore, that 
this second equilibration process deserves further study. Though the agreement in Case A 
(within a factor of three) has been considered good because of the uncertainties inherent at 
these low levels of radiation, this case also may benefit from further scrutiny of this process. 

Conclusions 

A study was conducted to evaluate the influence of several energy exchange models 
used in nonequilibrium computational fluid dynamics. Flowheld solutions obtained from 
LAURA with various forms of the exchange models were used as input to LORAN, a recently 
developed nonequilibrium radiation model. The results were evaluated by comparison to the 
Project FIRE II flight data taken at 1631, 1634, and 1637.5 sec in the early part of the 
trajectory. 

Models for the rate controlling temperature for dissociation, 2h, vibrational-translational 
energy exchange cross section, oy, and energy exchange in dissociation, A E v , were considered. 
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Formulations that increased 2h, reduced a v , or increased A E v resulted in reduced radiation. 
The greatest effect was seen at 1631 sec when nonequilibrium was most extreme. Radiative 
heating predictions from the different formulations varied more than 100 percent for this 
case. The lower radiation predictions generally fell closer to the flight measurements. 

In addition to examining the total radiation, predicted radiation spectra were compared 
with the flight data. These indicated that the close agreement in total radiation was fortu- 
itous in some cases because there was excessive energy in the molecular bands. At 1637.5 sec, 
a thermal equilibrium result provided better spectral agreement. At 1634 sec, an alternate 
chemical kinetics set improved the spectral comparison. These results suggest that chemical 
kinetic effects deserve further attention. 

This study has demonstrated that the examination of radiative heating measurements 
from flight experiments can provide a means of evaluating models used in computational 
fluid dynamics (CFD) for nonequilibrium flowfields. Further study of the FIRE radiometer 
data and results of future flight experiments can be expected to provide important tests of 
such models leading to improvements in CFD solutions. 
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Figure Captions 


1. Temperature Profiles for FIRE II - Case C 

2. Effect of Td Models on Molecular Dissociation - Case A 

3. Effect of Td Models on Temperature Profiles - Case A 

4. Effect of Td Models on Radiative Emission Profiles - Case A 

5. Effect of (T v Models on Temperature Profiles - Case A 

6. Effect of (T v Models on Radiative Emission Profiles - Case A 

7. Effect of Molecular Dissociation Models on Temperature Profiles - Case B 

8. Effect of Molecular Dissociation Models on Radiative Emission Profiles - Case B 

9. : a) Predicted Radiation Spectrum - Case A 

b) Measured Radiation Spectrum - 1631.3 sec 

10. : a) Predicted Radiation Spectra - Case B 

b) Measured Radiation Spectrum - 1634.43 sec 

11. : a) Predicted Radiation Spectra - Case C 

b) Measured Radiation Spectrum - 1636.43 sec 

12. Effect of Chemical Kinetics Model on Temperature Profiles - Case B 
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Table 3: Effect of Energy Exchange Models in Case A 


Run 

Flux, W/cni 2 (0.31-5.4 

eV) 

Identifier 

Atomic 

Mol. 

Atomic 

Total 


Cont. 

Band 

Line 


flight 




16 ± 20% 

vllO 

.IE-4 

.54 

.4E-3 

.54 

v!55 

.7E-5 

2.0 

.5E-2 

2.0 

< 

I— 1 

-a 

CO 

.IE-4 

1.1 

.9E-3 

1.1 

v273 

.IE-4 

3.2 

.2E-1 

3.2 

v373 

.IE-4 

.41 

.5E-3 

.41 

dk73 

.2E-3 

1.7 

.5E-1 

1.7 

v!732evs 

.IE-4 

.72 

.5E-3 

.73 



Figure 1: Temperature Profiles for FIRE II - Case G 
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cm 


Table 4: Effect of Energy Exchange Models in Case B 


Run 

Flux, W/cni 2 (0.31-5.4 

eV) 

Identifier 

Atomic 

Mol. 

Atomic 

Total 


Cont. 

Band 

Line 


flight 



8.2 ± 20% 

vllO 

.2E-1 

17.7 

2.7 

20.3 

v!55 

.4E-1 

19.1 

3.8 

22.9 

< 

I — 1 

-a 

oo 

.3E-1 

17.7 

3.3 

21.1 

v273 

.7E-1 

22.5 

5.2 

27.8 

v373 

.IE-1 

17.6 

2.2 

19.9 

dk73 

.3 

3.1 

7.9 

11.4 

v!732evs 

.1 

12.1 

6.0 

18.3 



y, m 


Figure 2: Effect of Tj Models on Molecular Dissociation - Case A 
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Table 5: Effect of Energy Exchange Models in Case C 


Run Flux, W/cm 2 (0.31-5.4 eV) 


Identifier 

Atomic 

Cont. 

Mol. 

Band 

Atomic 

Line 

Total 

flight 




81.7 ±20% 

vllO 

3.8 

33.5 

45.3 

82.7 

v!55 

3.9 

38.9 

46.2 

89.0 

< 

I — 1 

-I 

CO 

3.9 

35.4 

45.7 

84.9 

v273 

4.1 

38.8 

49.2 

92.1 

v373 

3.8 

33.6 

44.5 

81.9 

dkeq 

5.2 

18.8 

99.9 

123.9 

vleq 

6.5 

5.8 

75.6 

87.8 

vl732evs 

4.6 

20.8 

52.0 

77.4 

v!733dis 

4.6 

15.5 

52.6 

72.8 
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60000 r 



Figure 3: Effect of Tj Models on Temperature Profiles - Case A 



Figure 4: Effect of Tj Models on Radiative Emission Profiles - Case A 
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Figure. 7: Effect of Molecular Dissociation Models on Temperature Profiles - Case B 



y, m 


Figure 8: Effect of Molecular Dissociation Models on Radiative Emission Profiles - Case B 
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X, microns 


Figure 9: a) Predicted Radiation Spectra - Case A 


Figure 9: b) Measured Radiation Spectrum - 1631.3 sec 
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Figure 10: a) Predicted Radiation Spectra - Case B 


Figure 10: b) Measured Radiation Spectrum - 1634.43 sec 


0.60 
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Figure 11: a) Predicted Radiation Spectra - Case G 


Figure 11: b) Measured Radiation Spectrum - 1636.43 sec 
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Figure 12: Effect of Chemical Kinetics Model on Temperature Profiles - Case B 
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